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Abstract 

The system of second-order equations of mean-field expansion is considered for Nambu- 
Jona-Lasinio model with chiral symmetry of SU(2)-group. The system includes equations for 
four-particle and three-particle functions and also equations for next-to-leading two-particle 
function and next-to-next-to-leading quark propagator. Exact solutions for four-particle and 
three-particle equations are obtained. The solution for four-particle function is a disconnected 
combination of the leading-order two-particle functions. The connected part of three-particle 
function includes two-meson and three-meson contributions. The solution of three-particle equa- 
tion permits to close the equation for ncxt-to-lcading two-particle function and to calculate a 
correction to pion-decay constant This correction increases from 14% to 28% when values 
of quark condensate varies from — (210MeV) 3 to -(250AfeV) 3 . 

Introduction and Summary 

Nambu-Jona-Lasinio (NJL) model with quark content is one of the most successful effective models 
of quantum chromodynamics of light hadrons (for review see [l]-[3]). A number of physical appli- 
cations of NJL model is connected with multi-quark functions, which are the subject of present 
work. These multi-quark functions arise in higher orders of the mean-field expansion (MFE) for 
NJL model. To formulate MFE we have used an iteration scheme of solution of Schwinger-Dyson 
equation with fermion bilocal source, which has been developed in works [4]. We have considered 
equations for Green functions of NJL model up to the second order of MFE. The leading approxi- 
mation and the first order of MFE maintain equations for the quark propagator and the two-particle 
function and also the next-to-leading (NLO) correction to the quark propagator. The second order 
of MFE includes the equations for four-particle and three-particle functions and also equations for 
the NLO two-particle function and the next-to-next-to-leading (NNLO) quark propagator. 

We have found solutions for the four-particle and three-particle equations. The solution for 
four-particle function is a disconnected combination of the leading-order two-particle functions, 
consequently, the physical effects, which connected with four-particle function (i.e., pion-pion scat- 
tering), are suppressed in this order of MFE. 

The solution of four-particle equation gives us a possibility to close the equation for three- 
particle function. The solution of three-particle equation contains both disconnected and con- 
nected parts. The connected part of three-particle function includes two-meson and three-meson 
contributions. 

The solution of three-particle equation permits to close the equation for NLO two-particle 
function and to calculate a correction to pion-decay constant f n . This correction varies from 14% 
to 28% when values of quark condensate varying from -(2WMeV) 3 to -(250M e y) 3 . 

The correction to pion-decay constant has been calculated early in work [5] in the framework of 
improved l/n c expansion by using the partially bosonized version of NJL model (see also [6]). Our 
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result is in accordance with results of these work, though the method of the calculation is different. 
This fact demonstrates a self-consistency of NJL model concerning to high-order corrections. 



1 Mean-field expansion in bilocal-source formalism 

We consider NJL model with chiral symmetry SUy(2) x SUa(2). The model Lagrangian is 



C = ipidijj + — 



(1) 



Here tp is the quark field with n c colours, g is the coupling constant of m~ 2 dimension, r are Pauli 
matrices. 

The mean-field expansion in the bilocal-source formalism for the model can be constructed with 
the method of work [4] (see also [7] ) . 

Generating functional G of Green functions is the functional integral 



G{rj) = j D(ip, ip) exp ij J dxC — J dxdyifj(y)rj(y, x)ip(x)^ . 



(2) 



Here rj(y, x) is the bilocal source of quark field, n-th derivative of G over source rj is n-particle 
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As a consequence of translational invariance of the functional-integration measure in (2) we 
have the functional-derivative Schwinger-Dyson equation: 



8{x - y)G + id x 



5G 



5rj(y,x) 



+ ig\ 



5 



tr 



5G 



5r](y,x) \-5r](x,x)- 
5G 



75 T- 



5rj(y,x[ 



tr 



75T- 



5G 



5r/(x, x) 



(3) 



dx\rj{x, x\) 



5v(y,xi)' 

As a leading approximation for the mean-field expansion we consider equation (3) without right- 
hand-side. A solution of the leading-order equation is the functional 



G(°) = exp {Tr [s * 77] } 



(4) 



(Here and below Tr denotes the operator trace, and * denotes the operator multiplication.) Function 
S is a solution of the equation 



S(x) + idS(x) + igS(x) tr [S(0)\ = 0. 



(5) 



This equation is in essence the coordinate form of the well-known gap equation for the quark 
propagator. In the leading approximation the unique connected Green function is quark propagator 

S[ 0) =S = (m-p)-\ 

Here m is a quark dynamical mass. 

The leading approximation (4) generates the linear iteration scheme: 

G = G(°) + + • • • + G (n) + • • • , 

where functional G^ is a solution of the iteration-scheme equation 
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The general solution of equation (6) is the functional 

Q(n) _ p(n)Q{0) 

where is a polynomial of 2n-th degree on source 77. 



2 First-step equations 

The functional of the first step of the iteration scheme is 



-Tr 



S 2 * 



+ Tr 



Here S 2 is the two-particle function, is the NLO quark propagator. Taking into account the 
leading-order equations we have the system of equations for S 2 and . 
The equation for S 2 has the form 



So 
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Here tr u denotes the trace, which includes the upper line of function S 2 . 

A solution of equation (7) is well-known. Usually such equation is solved in momentum space. 
We give the solution of this equation in coordinate space since a similar method will be used for 
solving the equation for the three-particle function. 

Let us go to the amputated function 



and define the connected part 

x y 



F 2 = S- 1 * S- 1 * So * S- 1 * S" 1 
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Then for the function F% we obtain the equation 

x y 



n y x , *, j = -igS(x - y)S(x - y')S(x' - y){\ ■ 1 - 75 r- 75 r}+ 

+igS(x-y) J dzidyij tr u [S(x - xi)F% ^ X J ^ ^S(yi-y)]- 

-KTtv u [ l5 TS(x - xx)F 2 c (^ X x ) y y )^j S( V1 - y)}}. 

Iterations of this equation reproduce the discrete algebraic structure (colour, isotopic and Lorentz) 
of the inhomogeneous term. This circumstance permits to write the general form of the solution 
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Here A a is the scalar (sigma-meson) amplitude, A w is the pseudoscalar (pion) amplitude. These 
amplitudes are solutions of equations 

A a (x) = —ig5{x) + ig J dx\ tr [S(x — x±)S(xi — x)]A a (xi), 



In momentum space 



where 
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(10) 



(2vr) 4 (m 2 - {p + q) 2 ){m 2 - q 2 )' 

The integral in equation (10) is divergent and should be dealt with by some regularization. 

The equation for the NLO propagator is a system of simple algebraic equations. Introducing 
the first-order mass operator with formula = S^ 1 we obtain for the following 

expression 

^ 1 \x)=ig6(x)tv[S^(0)] + S(x)A (7 (x) + 3 l5 S(x) l5 A 7r (x). (11) 

3 Second-step equations. Solutions for four-particle and three- 
particle functions 

The second-step generating functional is 

G(2) = {hH s * * ^) + * + * ^) + H s(2) * ^) } G(0) ' 

i.e., second-step equations define four-particle function S4, three-particle function 53, and also NLO 
two-particle function and NNLO propagator S^ 2 \ For these four functions we have a system 
of four integral equations. All these equations (and all equations of following steps of the iteration 
scheme) possess the common structure, which is similar to the structure of equation (7): 
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Difference becomes apparent in the structure of inhomogeneous terms. The inhomogeneous term 
in the equation for S4 is 
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where £2 is defined in preceding section, i.e. the equation for S 4 does not include another three 
second-step functions. The inhomogeneous term in the equation for three-particle functions^ in- 
cludes function S 4 . This inhomogeneous term has the form 



-S(x 1 -y2)S{x 2 -yi)S^\x 3 -y 3 ) - S(xi - y 3 )S(x 3 - 2/i)S (1) (x 2 -2/2)- (14) 
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Similarly, the inhomogeneous term in equation for includes function S3, and the inhomo- 
geneous term of the equation for includes function . 

Due to such structure of the system its solution should be started by equation for four-particle 
function S4, then should be solved the equation for three-particle function S3, etc. 

The equation for four-particle function S4 with the inhomogeneous term (13) has the simple 
solution (see also [8]) 
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(15) 



This solution is disconnected, and it means the absence of physical effects due to four-particle 
functions in the given order of mean-field expansion. Particularly, the pion-pion scattering is 
absent in the given order and will appear in the next order. 

Solution (15) of the equation for four-particle function gives us the closed equation of type (12) 
for three-particle function S3. The solution of this equation with inhomogeneous term (14) can be 
obtained likewise to solving of equation (7) for the two-particle function. We go to the amputated 
function 

F 3 = S- 1 * S" 1 * S- 1 * S 3 * S~ x * S- 1 * S- 1 

and then separate the functions with the algebraic structures, which are reproduced after iterations. 
Connected part of the amputated three-particle function possesses two-meson and three-meson 
contributions. The explicit form of this function in momentum space see in Appendix. 

4 Two-particle function and pion-decay constant 

The solution of the equation for three-particle function permits to close the equation for NLO 
two-particle function ■ This equation differs from equation (7) only by inhomogeneous term. 
The inhomogeneous term in equation for s!p is 
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The two-particle function of the second step enables to calculate a correction to pion-decay constant 
/„-, which is one of the basic strong-interaction parameters. The pion-decay constant is defined by 
the relation 

iU5 w P^=<0\J%\P,b' >, 

where \P, b' > is the pion state with momentum P and isospin b, and J^ 5 = ^7^75 -yV' is the axial 

current. If the two-particle function has a pole term Sf ' 6 , which corresponds to the pion, then, 
taking into account these definitions, this pole term is connected with the pion-decay constant by 
relation 

~S(P ~ P')fl = * J dxdx'e^- p '^ t. u4 {wl ■ W\SI° 1 * ( * * ^] (16) 

Here tv u ^ denotes the traces over up and down lines of function S^ ' 6 . 

In the leading order we obtain from (16) the well-known expression for the pion-decay constant 
(see, for example, [1]): 

(/ (0))2 = _4in c m 2 / (0), (17) 

where Jo is defined by formula (10). For a regularization with four-dimensional cutoff 
Io(0) = j^jz log A2 +r 2 - K if m i , where A is the cutoff parameter. 

In the next-to-leading order formula (16) defines correction (/i^) 2 to expression (17). Surely, 
to calculate this correction it is no need in a complete solution of the NLO two-particle equation. 
In correspondence with (16) it is quite enough to calculate the pion-pole part only. 

The results of calculation of ratio rr = u T n / in the scheme with four-dimensional cutoff are 

shown in Table for three sets of model parameters. These sets are correspond to different values of 
quark condensate c = (i < ^pip >) 1 / 3 at physical value of pion-decay constant = 93 MeV. 



c (GeV) 


m (GeV) 


A (GeV) 


k = 3 5 A 2 /2tt 2 


Tf = (/i 1} ) 2 /(/^) 2 


-0.21 


0.42 


0.73 


1.87 


0.14 


-0.23 


0.28 


0.87 


1.33 


0.20 


-0.25 


0.24 


1.03 


1.19 


0.28 



Table. Values of ratio 77 = (/i^) 2 / (.fi°^) 2 in scheme with four-dimensional cutoff for different sets of values 
of chiral quark condensate c, quark mass m, regularization parameter A and dimensionless coupling k. 



Apparently the correction to value of the pion-decay constant increases from 14% to 28% when 
the absolute value of the quark condensate increases from (0.21 GeV) 3 to (0.25 GeV) 3 . These 
values of the correction to pion-decay constant are not far from results of calculations in works 
[5]-[6], though the method of calculations is different. This principal correspondence seems to be a 
reflection of some stability of NJL model with regard to quantum corrections. 
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Appendix 

The connected part of the three-particle funcion in the momentum space is 

pc ptwo— meson , pthree— meson 

3 3 3 

ptwo-meson = + p/j + + p//^ . A ^p/) . Aff (p")_ 

-A CT (P) • (S(p' + P") + 5(p' + P)) • A CT (P")- 

-A ff (P) • A CT (P') • (S( P " + P) + S{p" + P'))+ 

+( l5 S(p + P')75 + j 5 S(p + P")tb) • MP 7 ) • A.(P")+ 

+A CT (P) • ( 75 S(p' + P") + S(p' + P) 75 )r- A W (P")+ 

+A CT (P) • A^P') • r(S(p" + P) 75 + 7 5 S(p" + P'))+ 

+ (S( P + P')75 + 75^(P + P"))t ■ K(P') ■ A W (P")+ 

+A 7T (P) ■ ( 75 S(p' + P")75 + lsS{p' + P) 75 ) • A W (P")+ 

+A W (P) • A CT (P) • r( 75 S(p" + P) + 5(p" + P') 75 )+ 

+ ( 75< S(p + P') + S(p + P") 75 )r- A 7r (P') • K(P")+ 

+A n (P) ■ r{S{p' + P")75 + ^S(p' + P)) ■ A CT (P")+ 

+A 7r (P) • A 7r (P / ) • (75^" + P) 75 + 75^" + P')75)+ 

+i( 75 5(p + P')75 - 7 5 S(P + P")75)T- [A w (P) x A 7r (P // )]+ 

+iA 7r (P) • [( 75 5(p' + P")75 - 7 5 S(p' + P) 75 )r x A 7r (P // )]+ 

+iA w (P) • [A^(P') x T( 75 S(p ,/ + P) 75 - 75^(p" + P')75)]- 

Pf ree - meso " = 2n c A CT (P) • A CT (P') • A a (P")[A sss (P; P + P') + A SSS (P; P + P")]- 

-2n c k a {P) ■ A 7r (P / ) • A^(P")[A spp (P; P + P) + A spp (P; P + P")]- 

-2n c A n (P) ■ A n (P>) ■ A a (P")[A spp (P"; P> + P") + A spp (P"; P + P")]- 

-2n c A n (P) ■ A a (P') • A 7r (P ,/ )[A spp (P';P / + P") + A spp (P';P + P')]. 

Here p = p Xl , P = p Xl + p^ , p' = p X2 , P' = Pa; 2 + P V2 , p" = Px 3 , P" = Px s + Py 3 , where p Xi , p Vx are 
quark momenta and P + P' + P" = 0. We also use the notations A CT = 1 o A a for the scalar-singlet 
function, A n = 75T0 for the pseudoscalar-isovector function (see (9)) and 

A SSS (P;Q) = J dqtr[S(P + q)S(Q + q)S(q)}, A spp (P;Q) = J dqtr [S(P + q) 75 S(Q + q) 75 S(q)} 

for quark triangles. (In last formulae the traces are taken over Lorentz spinor indices only.) 
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